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ABSTRACT 

Self-similar shock solutions in spherically symmetric polytropic gas flows are con- 
structed and analyzed in contexts of proto-star formation processes. Among other 
possible solutions, we model a similarity shock across the sonic critical curve with 
an inner free-fall core collapse and a simultaneous outer expansion of the extended 
envelope; the separation or stagnation surface between these two flow zones travels 
outwards in a self-similar manner at a variable speed. One readily obtains accretion 
shock solutions. Semi-complete self-similar solutions across the sonic critical curve ei- 
ther once or twice without shocks can also be constructed. Features of star formation 
clouds of our polytropic model include the mass density scaling in the outer flow zone 
p oc r~ 2 /( 2 ~ 7 ), the temperature scalings of the inner flow zone T oc r -3 ^ 7-1 ^ 2 and of 
the outer flow zone T oc r^ 7-1 ^ 2-7 ', and the variable central mass accretion rate 
M = k^H^-^mo/G where 7 is the polytropic index, A: is a constant, too is the core 
mass, and G is the gravitational constant. Spectral line profiles characteristic of the 
'envelope expansion with core collapse' (EECC) shock solutions are expected. Ran- 
dom magnetic field permeated in a partially ionized cloud can be incorporated into 
this theoretical polytropic model framework. We discuss briefly our results in context 
of the oft-observed starless B335 cloud system as an example. 

Key words: ISM: H II regions — hydrodynamics — ISM: clouds — shock waves — 
stars: formation — stars: winds, outflows 



1 INTRODUCTION 

In astrophysical and cosmological contexts, the dynamic 
evolution of a spherical gas flow has been studied both 
numerically and analytically since late 1960s ( L arson T9691 
IPenston 19 69). Sufficiently far away from initial and bound- 
ary conditions, such flows may gradually evolve into a 
self-similar phase (Sedov 1959). Shu (1977) constructed 
the 'expansion wave collapse solution' (EWCS) in a self- 
gravitating isothermal gas cloud and advocated the so- 
called inside out collapse scenario for forming low-mass 
stars IShu 19771 |Shu, Adams fe Lizano 1987| . Such a sim- 
ilarity isothermal cloud collapse process involves a diverg- 
ing mass density p oc r~ 3 ^ 2 and a diverging radial in- 
fall speed v oc — r _1 ' a near the collapsing core yet with 
a p oc r~ 2 far away in a static envelope characterized 



by the outer part of a singular isothermal sphere (SIS). 
The boundary of the infall region expands at a constant 
sound speed. For a more general polytropic gas, similar- 
ity solutions have also been derived and analyzed from 
various perspectives l |Cheng 1 978 Goldrcich & Webe ~980l 
lYahil 19831 ISuto fc Silk l~988l |McLaughlin fc Pudritz 1997| 
|Fatuzzo, Adams fc Myers 2004| |Lou fc Wang 2006) . With 
effects of unknown energetic processes relegated to 
the polytropic index 7, these polytropic solutions, 
among others, become more flexible and versatile in 
modelling core collapsing processes in a molecular 
cloud. Shocks can be constructed in an isothermal 
flow (ITsai fc Hsu 19951 IShu et al. 20021 IShen fc Lou"2004l 
IBian fc Lou 20051 |YuT Lou, Bian fc Wu 2006| >. These simi- 
larity shock solutions contain important realistic features 
to model large-scale flow dynamics around star-forming re- 
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gions. While there are two- or three-dimensional numeri- 
cal simulations to model sub-structures in star-forming re- 
gions (True-lov e et al. 1998L we explore in this paper the full 
extent of possible one-dimensional features under spherical 
symmetry that are extendable to two- or three-dimensional 
flow modelling in principle. 

Recently, Lou & Shen (2004) obtained a novel class 
of semi-complete isothermal flows, including those referred 
to as 'envelope expansions with core collapses' (EECC) 
that may or may not involve the sonic critical line. 
Such solutions can model simultaneous inner collapse and 
outer flow expansion with or without shocks and are ap- 
plied to star-formation regions such as B335 cloud system 
llZhou et al. 19931 ISTien fe Lou 2004t . In particular, Shen & 
Lou (2004) emphasized that a cloud envelope, with a cen- 
tral collapsing core to form a star or a star cluster, may 
either expand or contract depending on initial environmen- 
tal conditions. Fatuzzo et al. (2004) considered a polytropic 
gas collapse under self-gravity in reference to observationally 
inferred far-away inflows. They adopted different polytropic 
indices for equations of state for an initial static cloud and 
a cloud in dynamic evolution in a more general manner. 
Central mass accretion rates are also examined as compared 
with observations (Fatuzzo et al. 2004). In this paper, we 
study polytropic gas flow solutions across the sonic critical 
curve either with or without shocks, which was not con- 
sidered by Fatuzzo et al. (2004). Moreover, semi-complete 
EECC solutions crossing the sonic critical curve twice in a 
polytropic gas are also constructed as an extension of the 
isothermal result of Lou & Shen (2004). Both shock and 
EECC features are expected in clouds with collapsing cores 
surrounding proto-star forming regions. This theoretical de- 
velopment, together with shock modelling (Bian & Lou 2005; 
Yu & Lou 2005; Yu et al. 2006) and adaptations, may also 
bear physical implications for the dynamic evolution phase 
with a timescale of ~ 10 yrs linking AGBs and PNe (e.g., 
Balick & Frank 2002 and extensive references therein) 

Observations of star-forming regions have been ac- 
cumulating over years. Various optical, infrared, sub- 
millimeter and more recent X-ray observations reveal struc- 
tures on different scales of hundreds of star-forming re- 
gions, a few of which are well observed l)Zhou et al. 19931 
IShirley fc Evans II 200*2] ISwift et al. 200511 and become 
valuable sources for checking and testing theoretical model 
development. By simulating characteristic spectral line pro- 
files given an underlying model of dynamic or magnetohy- 
drodynamic (MHD) flows, one can infer information of a 
star-forming region by estimating the density and flow pro- 
files, velocity dispersion, temperature scaling and central 
mass accretion rate and so forth. 

With these considerations in mind, we present in §2 our 
polytropic similarity flow model, including the polytropic 
EECC and shock solutions. Several aspects of probing star- 
forming processes are discussed in §3. Further discussion of 
our model application is presented in §4. Several details of 
mathematical derivations are summarized in Appendices A, 
B and C for the convenience of reference. 



2 SELF-SIMILAR POLYTROPIC 
GAS FLOWS 

For the dynamical evolution of spherically symmetric gas 
flows under self- gravity and thermal pressure, we use the 
basic nonlinear fluid equations in spherical polar coordinates 
(r, 9, (f>), namely 

d 
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Here mass density p, thermal gas pressure p and radial flow 
speed u are all functions of both radius r and time t, and 
equations Q and @ represent mass and momentum con- 
servations, respectively; M is the total enclosed mass inside 
r at a given time t and the first of equation ^ is another 
form of mass conservation; the Poisson equation relating gas 
mass density and the gravitational potential is automatically 
satisfied and G is the universal gravitational constant. 
The generalized polytropic relation takes the form of 

P = K(t)p-> , (4) 

where 7 is the polytropic index and the coefficient K(t) 
can vary with time t. The isothermal case corresponds to 
7=1 and a constant K(t) = k. Equation of state (|1J 
contains information of some unknown energetic processes. 
A time-dependent K(t) leads to different solution families 
for 1 < 7 < 2, while for p = hip 1 , the special cases of 
7 = 4/3 and 5/3 should be handled separately < |Cheng 1978| 
IGoldreich fc Weber 1980l|Lou fc Wang 2006) . 

We introduce the self-similarity transformation in the 
following form (ISuto fc Silk 19881 |Lou fc Wang 2006) 

x = r/A(t) , M(r,t) = B(t)m(x) , u(r,t) = C(t)v(x) , 



p{r, t) = D(t)a(x) , p(r, t) = E{t)(3{x) 



(•») 



where the forms of functions A(t), B(t), C(t), D(t) and E(t) 
should be properly chosen. Depending on x only, functions 
a(x), v(x), m(x) and /3(x) are the reduced density, radial 
velocity, enclosed mass, and thermal pressure, respectively. 
The scale- free independent variable x is defined following the 
power-law time dependence assumption A(t) oc t n with n 
being a constant index; this self-similar independent variable 
is thus 



x = r/{k 1/2 t") , 



(6) 



where k is a constant related to the polytropic sound speed 
c ~ fc 1//2 t n_1 . Here, the scale-free and self-similar indepen- 
dent variable x is made dimensionless and defined by x = 
r/(ct), and the time variation of the polytropic sound speed c 
is related to the power-law time dependence x = r/(k 1 ^ 2 t"). 

By taking into account of the following scaling relation- 
ships u ~ c, GM/r ~ u 2 , M ~ pr 3 , p ~ c 2 p and the 
equation of state we can consistently cast the variables 
into the following form (see Appendix A for a more detailed 
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discussion) 

p = a{x)/{4-KGt 2 ) , M = k 3/2 t 3n - 2 m{x)/[(3n - 2)G] , 

u = k 1/2 t n ~ 1 v{x) , p = kt 2n ~ 4 a 1 {x)/{A-KG) , (7) 

where, in reference to equation of state Kit) = 

fc(47rG) 7_1 t 2(7+n " 2) is readily identified, and the polytropic 
sound speed defined by 

c = (d P /8p) 1/2 = {ky) 1/a t n - 1 a^- 1 '> /a (8) 

varies with t and r in general. By the definition of K(t), 
K(t) = k becomes a constant coefficient for n = 2 — 7 and 
this corresponds to a conventional polytropic gas. 

Introducing the self-similar transformation J7J to non- 
linear partial differential equations Q, @, ((SJ, and J1J, we 
readily obtain two coupled nonlinear self-similar ordinary 
differential equations (ODEs) (see Appendix B for a more 
detailed derivation), namely 

[(n — nx) 2 — 7a 7 " 1 ] (dv/dx) = (v — nx) 2 a/(3n — 2) 

+2(v — x)ja 1 ~ 1 / 'x — (v — nx)(n — i)v (9) 

and 

\(v — nx) 2 — 7a 7 " 1 ] \da/(adx)] = (n — l)v 

— (v — nx)a/(3n — 2) — 2(v — x)(v — nx)/x . (10) 

In this procedure of derivation, an important result is 
m(x) — ax 2 (nx~v) (see equation lB6L giving a clear physi- 
cal requirement of nx—v > for a positive enclosed mass. All 
the solutions seriously considered in this paper must satisfy 
this necessary requirement. By setting v — for all x > 0, we 
obtain a static solution from equations © and II lit , namely 

v = , a = [2 7 (3n - 2)/n 2 ] 1/(2 " 7) x ~^-i) > 

m = (2 7 ) 1/(2 " 7) (3n - 2 ) 1/(2 - 7) n- 7/(2 - 7 'x (4 - 37)/(2 - 7) (11) 

(Suto & Silk 1988; Fatuzzo et al. 2004), referred to as the 
singular polytropic sphere solution (Lou & Wang 2006). 

We consider asymptotic solutions of equations @ and 
11011 for x — > + , x — > +00 and x near the sonic critical curve 
in the next subsection §2.1. Semi-complete solutions of the 
EECC type and shock solutions are constructed numerically 
in subsections 82.2 and 52.3. 



2.1 Sonic Critical Curves and Asymptotic 
Behaviours of Self-Similar Solutions 

Physically, when the travel speed of disturbances relative 
to the local flow speed is equal to the local sound speed, 
a singularity or a critical point arises. In steady flows of 
Bondi (1952) accretion and Parker (1958) solar wind, hy- 
drodynamic solutions are obtained to go across the sonic 
critical point smoothly. Mathematically, we can construct 
the sonic critical curve by simultaneously setting the coef- 
ficients of dv/dx, of da/dx, and of the right-hand sides of 
equations and (HOP to vanish. The mathematical form of 



the sonic critical curve is then specified by 

/— (t-1)/2 

nx - v = V7 Q o 
nxo — «o _ „ (xo — v ) (nx — v ) 



(n - l)v + 



-a -2- 



= 



(12) 



(3n — 2) xo 

where the subscript in variables xo, vo and ao indicates 
the explicit association with the sonic critical curve. The first 
equation in equation 11211 is obtained from the coefficients of 
dv/dx and da/dx. The second equation in equation 1121 can 
be obtained either from the right-hand side of equation 1101 
or from the right-hand side of equation JUJ , together with the 
first equation in equation 11211 . Several sonic critical curves 
of different parameters are shown in Figure 1 as examples 
of illustration. 

To clarify solution behaviours near the sonic critical 
curve and to construct solutions across the sonic critical 
curve, we study asymptotic solutions (first-order Taylor ex- 
pansion) of the form x = xo + 8, v = vo + vi5 and 
a — ao + aiS, where 8 is a small displacement in x. By 
substituting these expressions into equations @ and IIUH . 
the leading terms a\ = da/dx and v\ = dv/dx satisfy the 
following two equations 

V7^' 3)/2 



Ql 



Vi 



■ 2vo/x , 



(7 + l)v 2 + [4(7 — l)vo/xo + n — 4ry + l] Vi 
+2(2 7 - l)vl/xl + 2[a /(3n - 2) - 2n - 4 7 + 4]v /x 

+(n - 2)/(3n - 2)a + (2n + 4 7 - 4) = , (13) 

respectively. Once parameters n and 7 are specfied and xo, 
vo and ao are determined from equation 1 12t for the sonic 
critical curve, vi = dv/dx can be solved from quadratic 
equation 1131 : a\ = da/dx can then be obtained from the 
first equation in equation 11 1 Ml accordingly. These eigenso- 
lutions in the vicinity of the sonic critical curve can be uti- 
lized to construct analytically smooth semi-complete solu- 
tions across the sonic critical curve without a shock. As the 
second equation in equation 113H is quadratic in V\, along 
certain segments of the critical curve, the eigensolutions may 
not exist. This situation does arise in constructing solutions 
shown in Figure 2. 

To leading orders, asymptotic similarity solutions of 
coupled nonlinear ODEs © and 1 1 (H for x — > +00 and 
x — > + are summarized below. In the limit of x — » +00, 
we have 

nA 



a = Ax 



-2/r, 



2-yA 



7-1 



„(n-2)/n 

(3n - 2) 

E (2-2 7 -n)/n + Ba .(n-l)/n . ^ 



(2n + 2 7 - 3)n 
in the limit of x — > + , the free-fall solution is 
(2m ) 1/2 (3n-2) (2m ) 1/2 

» = ^175-. a = 

and in the limit of x - 
L-P type solution is 

x -> + : v = 2x/3, 

(16) 

where A, B, mo and a* are relevant constants of integration. 



2 x 3/2 > m = m ; (15) 

■ + , the polytropic counterpart of 

a — a* , m = a* x 3 (n — 2/3) , 
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Figure 1. Sonic critical curves with 7 = 1, 1.1, 1.2 and 1.3 
in the case of n = 2 — 7 for a conventional polytropic gas. The 
sonic critical curves trace the locus of the sonic critical points 
u = nr/t — c. The sonic critical curve becomes a straight line only 
for 7 = 1 corresponding to the isothermal case. With 7 = 1.1, 1.2 
and 1.3, the sonic critical curves diverge for very small x. 



These asymptotic solutions are found by comparing the lead- 
ing terms of nonlinear ODEs © and HUH for x — > +00 and 
x — > + , respectively. For x — » + , two different solutions for 
the cloud with or without a central core mo are considered. 
If needed, higher-order expansion terms in these asymptotic 
solutions can be readily determined. 

These leading scalings for x — > + [i.e., inner solutions 
for free falls 11511 and expansions 1161 1 are independent of 7. 
For x — > +00, there is a constant speed solution for n = 1 
corresponding to the isothermal case (Whitworth & Sum- 
mers 1985; Lou & Shen 2004). For x — > +00 and n = 2 — 7, 
the gas flow Mach number M = u/c remains constant as 
shown in equation 117H below. 



M 



Bx (n ~ 1)/n + 0[x {n ~ 1)/n ] 



C 7 1 / 2 a (7-l)/2 7 l/2,4(7-l)/2. r (l-7)A 



B 



y/2 ,4(7-1)/: 



+ 0(x 



-l/r, 



(17) 



2.2 Polytropic EECC Similarity Solutions 

We numerically solve coupled nonlinear ODEs © and Hot 
using available asymptotic solutions as 'boundary condi- 
tions'. Eigensolutions near the sonic critical curve are de- 
termined by equation 11311 in order to cross the sonic critical 
curve smoothly. However, for certain segments of the sonic 
critical curve, such eigensolutions of vi and ai do not exist. 
It is still possible to construct shock solutions across these 
sonic critical line segments (sec H2.3I I, 

It is simple to construct solutions not crossing the sonic 
critical curve by simply integrating from a chosen bound- 
ary condition. For solutions crossing the sonic critical curve 
smoothly, matching procedure of Lou & Shen (2004) is 
adopted for constructing semi-complete smooth solutions. 
That is, for the same solution, integrated results from the 



inner and outer boundary conditions should match with each 
other at a chosen meeting point. 

Displayed in Figure are different types of flow solu- 
tions crossing or not acrossing the sonic critical curve. In 
this —v versus x presentation, the first quadrant with x > 
and — v > physically represents collapse or inflow solu- 
tions, while the fourth quadrant with x > and — v < 
physically represents expansion, wind or outflow solutions. 
The solid curve spanning these two quadrants is the sonic 
critical curve as determined from equation 1121 . Three types 
of flow solutions are collapse/inflow solutions, envelope ex- 
pansion with core collapse (EECC) solutions and L-P type 
solutions are shown together in Fig. |5] 

These different kinds of dynamic solutions originate 
from different initial states, carrying sensible physical inter- 
pretations. Those curves diverging at small x have free-fall 
features, and the mo values marked along these curves rep- 
resent their core masses. We observe that for a larger core 
mass mo, the outer part of the cloud tends to be an inflow. 
The solution curve marked with mo = 3.20 is connected to 
an outer initial state of B — —0.1, which represents an in- 
flow (see equation 1141 and note that in this 7 — 1.2, n — 1 
case, B is the leading term of t)+oo). The flow curve with 
mo = 1.77 and B — 0.1 shows a breeze for x > 17.5, while 
the flow curve with mo = 0.66 and B = 3.572 shows an 
outflow at large x. On the other hand, both the inflow and 
outflow initial states at large x may involve an inner state 
without a central core (e.g., solution curves with B — —1.797 
and B = 2.230 in Fig. [3). In short, we have the capability of 
constructing self-similar flow solutions passing through the 
sonic critical curve and making it possible to connect dif- 
ferent inner and outer initial states in all feasible manner. 
Our polytropic model framework can accommodate various 
conditions for the inner and outer parts of a star-forming 
cloud. 

Solutions spanning both the first and fourth quadrants 
represent the EECC feature. We see that both inner core 
collapse and outflow features exist in a semi-complete EECC 
similarity solution. This may be heuristically viewed as a 
combination of Bondi (1952) accretion and Parker (1958) 
solar wind such that the sonic critical curve are smoothly 
crossed twice (Lou & Shen 2004). We see also that a L-P 
type solution and an inner collapse solution with mo = 0.40 
do not cross the sonic critical curve because of the absence 
of eigensolutions near the sonic critical curve (see equation 
1131 and the relevant discussion). Nevertheless, this problem 
can be solved by introducing a shock as shown in the next 
subsection. We focus on astrophysical applications of various 
EECC solutions in S3. 



2.3 Various Similarity Shock Solutions 

Shocks are introduced to match solutions from a given outer 
boundary condition and a certain solution from inner crit- 
ical points or inner boundary conditions. For a polytropic 
gas, the shock jump conditions in the shock front reference 
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Y=1.2,n=1 




Figure 2. Three types of polytropic flow solutions for 7 = 1.2 
and n = 1 are constructed by specifying relevant parameters. 
The solid line is the sonic critical curve. Dashed curves in the 
first quadrant are collapse and inflow solutions. The incomplete 
dashed curve in the fourth quadrant which does not cross the sonic 
critical line is a counter part of the L-P solutions. The dash-dotted 
curves span two quadrants are EECC solutions. 



framework are (Landau & Lifsliitz ~1959t 

P2/P1 = Qi/fia = (7 + l)Ml/[(i - l)Ml + 2] , (18) 



pa/pi = 27X1/(7 + 1) - (7 - l)/(7 + 1) 



(19) 



T 2 /T 1 = [27M 1 -(7-l)][(7-l)M?+2]/[( 7 +l) 2 M 1 ] , (20) 

where M is the flow Mach number in our reference frame- 
work, u is the radial gas flow speed in the shock reference 
framework and subscripts 1 and 2 denote the pre-shock (up- 
stream) and post-shock (downstream) sides of the shock, 
respectively. The ideal gas law is used here and T represents 
the gas temperature. These shock conditions are derived 
from the mass conservation, the momentum conservation 
and the formal entropy conservation along the streamline. 

For gas flows with a constant polytropic index 7 for 
the equation of state, the sound speed c jumps across the 
shock front because of the temperature discontinuity there. 
That is, ca/d = (p 2 /p2) 1/2 /(pi/pi) 1/2 = (Ta/Ti)V a , where 
ci and C2 are pre-shock (upstream) and post-shock (down- 
stream) sound speeds, respectively. Our definition of the self- 
similar independent variable x is of the scaling x ~ f/{ct) 
(see equation |S| and the relevant discussion). Due to this 
jump in the sound speed c, independent variable x shifts 
in the self-similar framework accordingly across the shock 
front, namely 

xa/m = (Ti/T 2 ) 1/2 . (21) 

Therefore the same physical interface located at r at a given 
time t shows a separation across the shock front in our 
self-similar presentation because of this shock discontinuity 
shown in Figure [3] 



From equations Q and 1181 . we obtain the shock jump 
conditions in the shock reference framework in the reduced 
self-similar form 

01(7 + l)(«i — nn) 2 



Q.2 



(7 — l)( u i — nxi) 2 + 2-yaJ 



7-1 



V'Z 



nx2 + 



X-2 



353.(7+1) 



(7 — l)(vi — nxi) + 



7-1 -1 



(22) 



(23) 



vi — nxi 

(see Appendix C for details of derivation). Because of the 
entropy change across the shock front, the constant k is 
different on two sides of the shock front and the relation 



U/2 



X\jx% holds. We determine the post-shock 



quantities from the pre-shock physical conditions by apply- 
ing equations 122H and 12311 . During this process of shock 
construction, equations I20H and 12111 are used to determine 
the post-shock (downstream) independent variable xi. Thus 
these shock conditions serve as the condition to cross the 
sonic critical curve. Using the matching procedure of Shen 
& Lou (2004), we obtain different classes of shock solutions. 
That is, when integrated from an outer boundary condition 
to a certain place near the sonic critical curve, the jump con- 
ditions (1221 and i!2-Sfl are used to get the post-shock (down- 
stream) quantities. Simultaneously, another asymptotic so- 
lution from an inner boundary condition is integrated out- 
wards. When these two solutions match with each other in 
all quantities satisfying the shock conditions, a shock solu- 
tion is then found. 

Displayed in Figure are different types of shock so- 
lutions with polytropic parameters 7 = 1.2 and n = 0.8. 
In this presentation of —v versus x, the first quadrant with 
x > and — v > physically represents collapse or inflow 
solutions. While the fourth quadrant with x > and — v < 
physically represents expansion, wind or outflow solutions. 
The solid curve spanning these two quadrants is the sonic 
critical curve. Dash-dotted straight lines connect the pre- 
shock and post-shock sides across the shock. These straight 
connection lines are not vertical because the same shock 
front has different scalings for the independent variable x 
on two sides in the self-similar framework (see equation 
I21H . Three types of shock solutions are presented together, 
namely, shocks running into a static atmosphere, accretion 
shocks and wind shocks are shown in Figure 

We see that for the curve with the core mass param- 
eter mo = 1.51, different kinds of shocks (shocks running 
into a static atmosphere and wind shocks etc.) can be con- 
structed when the shock location or speed changes. There- 
fore the shock location is an additional parameter to join 
different inner (final) and outer (initial) conditions. For the 
same reason, shocks running into a static atmosphere (e.g., 
curves with xi = 2.06 and x\ — 4.4) are constructed in both 
clouds with and without a central core. 

From the two B = — 2 and B — 0.5 curves in Figure 
we note that an accretion shock can be constructed for a 
cloud either with or without a central core mass. A close look 
at the B = 0.5 curve shows that this is actually a 'breeze' 
shock solution which tends to be an outflow at a large x 
and the flow speed eventually approaches zero as x — > +00 
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Figure 3. Three types of shock flows are constructed for 7 = 1.2 
and n = 0.8 with other parameters specified for a conventional 
polytropic gas. The solid line spanning the first and fourth quad- 
rants is the sonic critical curve. Dash-dotted straight lines connect 
the pre-shock and post-shock sides (see equation I21l and the rele- 
vant discussion). The two dashed curves joining the — v = axis 
at x\ = 2.06 and x\ = 4.4 are self-similar shocks running into a 
static atmosphere (i.e., the outer part of an singular polytropic 
sphere). The two curves with outer initial states characterized 
by B = —2 and B = 0.5 show shock discontinuities in the first 
quadrant; these are examples of accretion shocks. In contrast, the 
curve with B = 1.23 describes a wind shock. 

(not shown in Fig.[3]because of scale restrictions). The pure 
coreless wind shock with B = 1.23 is a new addition to the 
shock solution family. 

These shock solutions show various features for the in- 
ner and outer parts as in Fig. [2] The introduction of shocks 
enriches the variety of the semi-complete solution family 
(because shock location or speed itself is an additional pa- 
rameter) and is qualitatively consistent with observations of 
star- forming regions (Ni sini et al. 19 991. 



3 APPLICATIONS TO STAR-FORMING 
CLOUD CORES 

There are several important observational aspects for an iso- 
lated star-forming region: the mass density profile, the tem- 
perature profile, the radial flow speed profile, the core mass 
accretion rate, magnetic field, and the angular momentum. 
In a spherically symmetric model, we consider and discuss 
the first five aspects except the last, and provide estimates 
of the star-forming cloud system B335. 

Before coming to specific astrophysical applications, we 
clarify the physical meaning of similarity solutions passing 
through the sonic critical curve. For recent observations of 
star-forming regions, the features of inner part and outer 
edge are well observed and may serve as initial states for 
different kinds of star-formation models. Similarity solu- 
tions passing through the sonic critical curve with or with- 
out shocks can join all kinds of different inner and outer 



physical states (see Figs.|5|and|SJl, while similarity solutions 
without crossing the sonic critical curve have much less free- 
dom (Suto & Silk 1988; Fatuzzo et al. 2004). For example, 
the EECC solutions (Lou & Shen 2004), most of which pass 
through the sonic critical curve, represent a new class of sim- 
ilarity solutions that are not constructed previously without 
encountering the sonic critical curve. Polytropic shock so- 
lutions also form a new family; these shocks connect var- 
ious inner and outer physical states and provide theoreti- 
cal basis for modeling shock phenomena observed recently 
(^Nisini et al. 1999^ . 

3.1 Mass Density and Temperature Profiles in a 
Cloud with a Free-Fall Collapsing Core 

By asymptotic solution 1151 at small x, the inner collapsing 
core in a free-fall state of a pre-stellar region has a mass 
density profile of p oc r~ 3 ^ 2 . This leading radial power-law 
scaling is also predicted by the isothermal model ( Sh u 19771 
ILou fc Shen 200411 and appears consistent with observations 
dZhou et al. 19901 [Andre et al. 2000l . In other words, this 
information of inner mass density scaling alone cannot dis- 
tinguish between polytropic and isothermal cloud models. 

For the outer part of a pre-stellar region (sufficiently far 
away from the collapsing core), the isothermal cloud model 
gives a mass density profile p oc r~ 2 , while the polytropic 
cloud model predicts p oc r _2// " (see equation I14H . Then 
for n — 2 — 7 corresponding to a constant K (£) = k in a 
conventional polytropic gas, the mass density profile for the 
outer part is 

p oc r-W-^ (24) 

So 7 > 1 means a density exponent smaller than —2 cor- 
responding to a steeper density profile, while 7 < 1 means 
a density exponent between —1 and —2 corresponding to a 
less steep density profile. Earlier observations gave a mass 
density profile with an exponent around —2 with consid- 
erable errors (IZhou et al. 19901 . Recent observations show 
that this mass density exponent for the very outer parts 
(at radii r > 15, OOOau) of starless cores is about —3 to 
—4 \ Andre et al. 2000)l . corresponding to a polytropic index 
7 ~ 1.3 — 1.5. This observational evidence appears sensible 
and encouraging for a polytropic model. 

Corresponding to the inner mass density profile p oc 
r _3//2 , the radial scaling of thermal temperature T for an 
ideal gas in the inner free-fall region would be 

T oc p/p oc r - 3(7 - 1)/2 . (25) 

This inner temperature T would be constant for an isother- 
mal cloud model with 7 = 1. The observed negative index 
of temperature scaling (Shir ley fc Evans II 2002| l calls for 
7 > 1 in the inner region, qualitatively similar to the 7 
value inferred for the outer mass density profile. These ra- 
dial scalings of cloud temperature are important features to 
constrain a polytropic cloud model by equations 1241 and 
d2~5l . 

For the reference of observational diagnostics, the gas 
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temperature sufficiently far away (i.e., large r) from the col- 
lapsing core is given either by a constant T for an isothermal 
cloud model or by 

T OC p/p OC r - 2 ("?-l)/(2-7) ( 26) 

for a conventional polytropic cloud model. Of course, the 
complete temperature profile in a conventional polytropic 
gas can be readily computed for observational comparison. 
The key message here is that independent determination of 
cloud temperature is extremely valuable for distinguishing 
models for cloud dynamics. 

3.2 Radial Gas Flow Speed Profiles 

Outflow features in star-forming regions are well docu- 
mented in recent years. In spite of the large-scale bipolar 
morphology of these outflows, we are mainly interested in 
core regions of these outflow sources with free-fall collapse 
or infall characteristics HAndre et al. 20001 IWu et al. 20041 
Swift laTal. 20051 that are grossly in accord with our shock 
scenario of EECC. 

Collapses and outflows were thought to happen in dif- 
ferent stages during a star formation process, and they were 
considered separately l)Shu 19771 [Yahil 1983). However, at 
least ~ 12% of identified outflows show characteristics of 
collapses or infalls llWu et al. 200411 . Note that infall char- 
acteristics are more difficult to identify than those of an 
outflow due to the complexity of observed spectral lines and 
in the analysis of line profiles. One would expect that more 
infall motions can be identified in outflow sources in the 
future. 

Our polytropic EECC shock solutions provide under- 
lying dynamic models for simulating spectral line profiles, 
during a certain epoch of a star formation cloud. The most 
expected line profile features are those from EECC solu- 
tions, that is, inflow and outflow within and outside a cer- 
tain radius r. Double peak spectral line features representing 
outflow will be observed away from the core, and combined 
double peak line features from both outflow and inflow will 
be seen around collapsing cores. 

Using an isothermal model for self-similar cloud dynam- 
ics (Lou & Shen 2004), Shen & Lou (2004) emphasized the 
notion that a molecular cloud with a free-fall collapsing core 
at the centre can have various possible forms of motion suf- 
ficiently far away from the collapsing core, depending on 
the evolutionary history of the cloud. Such far away radial 
flows of a cloud may be characterized by a wind, a breeze, 
an outer part of a static singular isothermal sphere (SIS), a 
contraction, or a inflow. Furthermore, all these possible ra- 
dial flow profiles can involve shocks across the sonic critical 
curve (Tsai & Hsu 1995; Shu et al. 2002; Shen & Lou 2004; 
Bian & Lou 2005). And this can give rise to a wide range 
for the central mass accretion rate to be discussed presently 
in the next subsection. All these important features in an 
isothermal cloud model can be further incorporated into a 
polytropic cloud model. 

In reference to extensive cloud observations, Fatuzzo et 



al. (2004) noted a ubiquitous feature of far away inflows 
and/or contractions in clouds with central free-fall collaps- 
ing cores. This differs from the static SIS envelope as initially 
introduced by Shu (1977) for the EWCS for star formation. 
The polytropic cloud model of Fatuzzo et al. (2004) also 
allows the polytropic indices to be different for an initially 
static cloud and for the subsequent evolution of a dynamic 
cloud. Nevertheless, Fatuzzo et al. (2004) did not consider 
self-similar polytropic flow solutions that cross the sonic crit- 
ical curve smoothly and that involve shocks. Our model anal- 
ysis here opens up these two important possibilities in a uni- 
fied theoretical framework. Moreover, we emphasize again 
that in a polytropic cloud model, a cloud with a central 
free-fall core collapse can be simultaneously characterized 
by a far-away polytropic wind (B / in equation I14|l or 
a far-away polytropic breeze (B = and A ^ in equa- 
tion with or without shocks. Gas flows with central core 
collapse and far-away outflow or breeze are common in star- 
forming regions <IWu et a l. 2004), indicating the capacity of 
our solutions passing through the sonic critical curve. 

It is highly desirable to synthesize spectral line profiles 
using radial flow speed and mass density profiles from these 
underlying self-similar dynamic models for developing radia- 
tive diagnostics and for extensive observational comparisons. 

3.3 Variable Central Mass Accretion Rate 

The core mass accretion rate from equations © , J7J and 1151 
is given by M = fc 3/2 i (3n_3) m /G. For n = correspond- 
ing to a constant K(t) — K of a conventional polytropic gas, 
the core mass accretion rate appears as 

M = fc 3/ Y 3 - 37) m /G (27) 

(see also equation 31 of Fatuzzo et al. 2004; the differ- 
ence in a numerical factor is related to a different definition 
of notations), which is constant in the isothermal case of 
7=1 ijShu 1977|l . A constant M is inconsistent with ob- 
servations within a formation timescale of ~ 10 5 — 10 6 yr 
\ Andre et al. 200(51 . The polytropic analysis does offer a 
possible leeway to model a variable M in time t. For a fixed 
mo, there are two parameters k and 7 to adjust. Here k is 
related to the sound speed c (see equation^ and 7 can be 
inferred observationally from the radial mass density pro- 
file as discussed in §3.1. For various shock solutions, the mo 
parameter can vary in a wide range. Although not easily 
done, it is crucial to identify the evolution phase of a dy- 
namical cloud system and estimate or guess t accordingly as 
discussed presently. 

For 7 < 1 and 7 > 1 in equation l|270 . the 
core mass accretion rate M increases and decreases with 
time t, respectively. Observations of star-forming systems 
llAndre et al. 2000 ) indicate that mass accretion rates gener- 
ally decrease in time t. This appears consistent qualitatively 
with 7 > 1 for most cloud systems with free-fall collapsing 
cores. This inference of 7 > 1 also agrees with the interpre- 
tation of steeper outer density profiles and the inner tem- 
perature power- law scaling with a negative index (see §3.1). 
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For shock flows, the core mass accretion rate varies with 
shock locations and thus shock speeds llShen fc Lou 20041 
IBian fc Lou 20051 , This is because the value of mo is differ- 
ent for different shock locations. 

The reason why 7 > f or 7 < f induce different M is as 
follows. For n — 2 — 7 in equation ©, a physical interface 
for a fixed x travels with increasing speed for 7 < f . When 
this interface is set between the inflow and outflow regions 
at v — 0, we have an increase of infall mass M oc pr 2 dr/dt. 
For a fixed x, the product pr 2 also increases with time t 
for 7 < f , then the increasing interface speed dr/dt leads to 
increasing mass infall rate. For 7 > f , the interface travels at 
a decreasing speed and the product pr 2 also decreases, thus 
the mass infall rate decreases. When we set the interface at 
the shock location, the conclusion is that 7 < f and 7 > f 
lead to acceleration and deceleration of shocks, respectively. 

3.4 Star-Forming Cloud System B335 

The cloud system B335 is a well-studied protostellar collapse 
candidate that has been observed through molecular lines 
l|Zhoii et al. f990l|Harvey et al. 2003|lEvans II et al. 2005) . 
Data fitting and analysis of H2CO and CS molecular transi- 
tion lines agrees well with the inside-out collapse model (Shu 
f977; Shu et al. f987), yet the angular resolution of obser- 
vations is still not good enough ijZhou et al. 1993^ . Recent 
observations of Institut de Radio Astronomic Millimetrique, 
Plateau de Bure Interferometer (IRAM PdBI) give several 
new results on this protostellar core ( |Harvey et al. 2003| > . 
First, a compact component with a circumstellar disc was 
identified. Secondly, the inner mass density profile is p oc r~ p 
with p = 1.55 ± 0.04 in a certain region r < Rq. Thirdly, 
another valuable information is the gas temperature profile 
of T oc r~ 0A in a certain region r < Rt- Here, we tentatively 
take R ~ Rt ~ 5000au. 

The inner mass density profile p oc r~ p with p — 1.55 ± 
0.04 agrees with a free-fall core collapse (p — 1.5). Both 
polytropic and isothermal models show this scaling result. 
The outer density profile (3, 500 — 25, OOOau) also shows a 
power law with an index p = 1.91 ±0.07. This result may be 
regarded as consistent with the isothermal model with p — 2, 
while the polytropic model can fit this result more accurately 
if we set 7 = 0.95 in the outer part. One difficulty for the 
isothermal model is that the shallow slope (p ~ 1) of density 
profile at the edge of collapse region is still not observed 
l |Harvey et al. 2003^ . But no shallow slope is necessary in 
the polytropic cloud model. 

The temperature profile is T oc 7- -0 ' 4 for r < Rt = 
5000au, that is T = 10K at r > R T and increase to T ~ 50K 
at r ~ lOOau {Harvey et al. 2003| > . As observations show 
that p oc r~ 1,55 near the core, we need to set 7 = 1.26 in 
the inner region to fit T oc r -0 ' 4 (see equation 1251 . While 
one might use the isothermal model and assume different 
temperatures in different cloud regions, the polytropic cloud 
model with a temperature variation appears more sensible. 

From the outer mass density profile and the inner tem- 
perature scaling of B335 system, we infer 7 = 0.95 for 



3500au < r < 25, OOOau and 7 = 1.26 for r < 5000au. 
This difference of 7 shows that 7 is not a constant at least 
for certain cloud systems with core collapses. Only when 
the outer mass density profile is steeper than p oc r~ 2 , can 
7 become a constant. 

Outflow and shock phenomena in B335 are also reported 
in recent years JNisini et al. 19991 1. While previous isother- 
mal solutions (Shu 1977) or polytropic solutions without 
crossing the sonic critical curve (Fatuzzo et al. 2004) cannot 
fit these observed conditions, these physical conditions can 
be related and constructed using an EECC shock curve pass- 
ing through the sonic critical curve. Thus, this cloud system 
may be modelled as EECC self-similar shock solutions in a 
polytropic gas. 

3.5 Magnetic Fields in Star-Forming Clouds 

When a cloud is partially ionized, the magnetic field per- 
meated in the cloud can exert an influence on the MHD 
evolution of the cloud system. Physically, magnetic field di- 
rectly interacts with the charged particles which in turn in- 
teract with neutral particles through sufficiently frequent 
collisions. If the MHD timescale of interest is much longer 
than collision timescales between electrically charged and 
neutral particles, we may ignore the effect of ambipolar dif- 
fusion and adopt a MHD approximation to describe large- 
scale flows of a magnetized cloud. Moreover, the presence 
of magnetic field also leads to useful radiative diagnostics. 
For example, the presence of relativistic electrons produced 
by shocks or proto-stellar activities and of magnetic field to- 
gether can give rise to cyclotron and synchrotron emissions. 

Depending on the evolutionary history, magnetic field 
in a collapsing cloud under self-gravity may possess a vari- 
ety of gross geometries (e.g., Chandrasekhar & Fermi 1953; 
Chandrasekhar 1954; Strittmatter 1966; Bisnovatyi-Kogan, 
Ruzmaikin, & Sunyaev 1973; Mestel 1986; Rees 1987; Lou 
1996). We outline below a possible scenario involving mag- 
netic field to pave the way for further investigations. It is per- 
ceived that a partially ionized gas cloud under self-gravity 
is permeated with a completely random magnetic field (e.g., 
Zel'dovich & Novikov 1971). For such a random magnetic 
field, we envision a random 'thread ball' scenario in a vast 
spatial volume of gas medium. A magnetic field line follows 
the 'thread' meandering within a thin spherical 'layer' in 
space in a random manner. Strictly speaking, there is always 
a random weak radial magnetic field component such that 
random magnetic field lines in adjacent 'layers' are actually 
connected throughout in space. By taking a large-scale en- 
semble average of such a magnetized gas cloud, we are then 
left with packed 'layers' of random magnetic field compo- 
nents dominantly transverse to the radial direction. Over 
large scales, the MHD collapse and flow in such a magne- 
tized cloud are approximated as quasi-spherically symmetric 
(Yu & Lou 2006; Yu, Lou, Bian & Wu 2006; Wang & Lou 
2006). 

In addition to the thermal gas pressure force and self- 
gravity in the cloud, the magnetic Lorentz force consisting 
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of magnetic pressure and tension forces will come into play 
for a sufficiently strong magnetic field. MHD generalizations 
of isothermal hydrodynamic self-similar solutions are readily 
available (Yu et al. 2006). For example, for the MHD free- 
fall asymptotic solution towards the core, the magnetic field 
strength diverges yet the self-gravity remains overwhelming 
over the thermal pressure force and the Lorentz force com- 
bined. In such an MHD core collapse process, a strongly 
magnetized environment can be created and sustained sur- 
rounding a proto-stellar core. Accretion and rebound MHD 
shocks can emerge (Bian & Lou 2005; Yu et al. 2006; Wang 
& Lou 2006). Near the very centre, various magnetic activi- 
ties, persistent or sporadic, are naturally expected to break 
the quasi-spherical symmetry as well as the self-similarity. If 
a cloud possesses an initial angular momentum, then such a 
core collapse under self-gravity will lead to a rapid gas rota- 
tion and a disc formation surrounding the proto-stellar core; 
in such a strongly magnetized environment, the coupling 
between the magnetic field and disc rotation can give rise 
to collimated outflows such as jets or winds (e.g., Lovelace 
1976; Blandford & Payne 1982; Shu et al. 1997; Ferrari 
1998). This type of MHD collapse processes may lead to 
intensely magnetized proto-stellar objects. 



4 DISCUSSION 

For star-forming cloud systems, the outer mass density pro- 
file, the inner and outer temperature profiles, and the core 
mass accretion rate are crucial to test or constrain a poly- 
tropic cloud model. Given a self-similar dynamic profile in a 
polytropic gas, spectral line profiles can be simulated. These 
spectral line profiles are very important to show different 
evolution epochs in forming proto-stars. In addition to form- 
ing proto-stars, the variable core mass accretion rate M by 
equation 1271 also bears physical implications for the dy- 
namic evolution phase with a timescale of ~ 10 3 yrs link- 
ing AGBs and PNe IIBalick fc Frank 200211 . Further obser- 
vations and analysis are needed especially for post-PNe of 
different ages. We also expect adaptations of our polytropic 
shock model to a certain epoch of supernova explosions (Lou 
& Wang 2006), with radiation and neutrino pressures also 
taken into account. 

Shock phenomena are common in forming proto-stars, 
formation and evolution of PNe, and supernova explosions; 
they are important for the gas heating, the production of 
high-energy photons, the synthesis of heavy elements, and 
the acceleration of cosmic ray particles. In these contexts as 
well as others, our polytropic shock model has wide appli- 
cability. Properties of polytropic shocks should be explored 
further to include radiation processes, magnetic field effects 
(Yu et al. 2006; Wang & Lou 2006) and atomic/nuclear pro- 
cesses associated with shocks. 

Radiative transfer processes will be incorporated into 
our model to make the underlying dynamic processes 'vis- 
ible'. Spectral line profiles of certain atoms and molecules 
can reveal profiles of mass density, flow velocity and gas 



temperature. The EECC profiles and shocks are the most 
important features. Certain spectral line features derived 
from various underlying dynamic solutions should be com- 
pared with observations. Recent and further high-resolution 
spectral IR and microwave telescopes will surely stimulate 
the development of dynamic study of star- formation, forma- 
tion/evolution of PNe, SNe and other relevant astrophysical 
processes. 
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APPENDIX A: DERIVATION OF 
SELF-SIMILAR TRANSFORMATION 

We here provide a heuristic derivation of self-similar trans- 
formation 10. In the end, it is the self consistency counts. 
The independent self-similar variable is [see eq. ijfjj] 

x = r/(k 1/2 t n ) (Al) 

and thus the polytropic sound speed c scales as ~ fc 1 ' 2 ^™ -1 . 
As the radial flow speed u and the polytropic sound speed 
c are expected to scale similarly, i.e., u ~ c, the radial flow 
speed takes the following form 

u = C(t)v(x) = fe 1/2 t n - 1 «(a:) . (A2) 

From the scaling relationships M oc pr 3 and GM/r 2 ~ 
u 2 /r, we have p related to u and r in the form of p ~ u 2 /Gr 2 . 
In reference to equations 1A1I and JA2I) . we then have the 



mass density p = v 2 (x)/(Gt 2 x 2 ). Therefore, we put the fac- 
tor 1/x 2 into the scale-free part and multiply a constant 
factor l/(47r) for notational convenience and arrive at 

p = D(t)a(x) = a(x)/(4ivGt 2 ) . (A3) 

There implies a relationship a(x) oc v 2 (x)/x 2 . 

From the scaling relationship M ~ pr 3 and the two 
self-similar expressions (E3 and l|A3jl . we have the enclosed 
mass M in the form of M = k 3/2 t 3n - 2 x 3 a(x) / AnG. We then 
merge the factor x 3 into the scale-free part and multiply a 
numerical factor 47r/(3n — 2) for notational convenience. It 
follows that the enclosed mass is in the form of 

M = B(t)m(x) = k 3/2 t 3n - 2 m(x)/[(3n - 2)G] . (A4) 

It is also apparent that the scaling m(x) ~ x 3 a(x) holds. 

The thermal gas pressure p has the scaling of p ~ c 2 p ~ 
kt 2n ~ 4 / (AtvG), where the form of the polytropic sound speed 
c oc fc 1 ^ 2 t n_1 and equation l|A3(l are used. For the scale-free 
part, we refer to the equation of state Q an d obtain the 
relation p = Kit^p 1 ~ oP(x). Therefore, the thermal gas 
pressure takes the following form 

p = E{t)/3(x) = fci 2n - 4 a 7 (a;)/(47rG) . 

It immediately follows that f3(x) ~ a 1 (x). 



(A5) 



APPENDIX B: DERIVATION OF 
SELF-SIMILAR ORDINARY DIFFERENTIAL 
EQUATIONS 

Here we provide the derivation of the two nonlinear coupled 
ODEs @ and 1101 using the self-similar transformation J7J|. 

By the chain rule, the partial differentials of a general 
function J-(r, t) = F(t)f(x) are: 

F(t) df(x) 



df(r, t) _ df{x) dx 



dr 



dx dr kMH n dx 



8F(r,t) _ dF(t) 4f(x) dx 



di 



dt 
dF(t) 



fix) 



dx dt 
nrF(t) df(x) 



(Bl) 



(B2) 



dt feVatn+i dx ' 

Referring to these partial differentials and the self-similar 
transformation £J, equations Q, and J5J can be readily 
reduced to 



, .da dv 2av 

(v — nx)— h a— — I 2a = 

dx dx x 



. . . . dv y _ 2 da 

(n — l)v + [V — nx)— h7« -; h 

dx 



dx (3n-2)x 2 



(B3) 



, (B4) 



(3n — 2)m = (nx — v)^^- , — — = (3n — 2)x 2 a , (B5) 
dx dx 

respectively. From equation HB5L it follows that 

m(x) = ax 2 (nx — v) , (B6) 

leading to a necessary physical requirement of nx — v > 
for a positive enclosed mass; this is referred to in Section 2 
of the main text. 
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Finally, from equations 1B3L 1B41 and 1B6L we derive 
the two coupled nonlinear ODEs in terms of v(x) and a(x) 

[(v — nx) 2 — 7a 7-1 ] (dv/dx) = (y — nx) 2 a/(3n — 2) 

+2(v — x)ya~ r ~ 1 /x — (v — nx)(n — l)v (B7) 

and 

\(v — nx) 2 — 7Q 7 ~ 1 ] [da/(ada;)] = (n — l)v 

-(v - nx)a/(3n - 2) - 2(v - x)(v - nx)/x . (B8) 



we obtain from equations (J7J and 1181 the density disconti- 
nuity in the self-similar framework as 



C*2 



an (7 + 



nxi ) 



(7 — l)(«i — nxi) 2 + 2-yaJ' 



(C7) 



Finally from equations 1181 , 1C3I — 1U6II , we derive the radial 
flow velocity discontinuity in the self-similar framework as 

2 7 g?-' 



nx2 + 



X-2 



Xl{j+1) 



(j~l)(vi-nx 1 ) + - 



Vi — nx\ 



(C8) 



APPENDIX C: DERIVATION OF SHOCK 
CONDITIONS IN THE SELF-SIMILAR 
FRAMEWORK 

We summarize here the derivation of shock conditions in the 
self-similar framework, namely, equations \2'2l and il'-MI) in 
the main text. 

In the reference framework of the shock front, the shock 
jump conditions for the mass density and flow velocity are 
shown in equation I18H . where subscripts 1 and 2 denote the 
pre-shock (upstream) and post-shock (downstream) sides, 
respectively. Note that flow velocity in the reference frame- 
work of the shock front has the form of Ui = Ui — u s and 
U2 — U2 — u s , where u s is the radial travel speed of the shock 
front. The shock front can be identified with x± in the pre- 
shock (upstream) side. Thus referring to equation ©, the 
travel speed of a shock front is 

u s = dr/dt — d{xik]/ 2 t n )/dt = nk^J 2 ^^ 1 x\ . (CI) 

This travel speed of the shock front can be also derived from 
the post-shock (downstream) side, namely 

u s = dr/dt — d[x2k\ t n ) f dt = nk^ 2 t n ~ 1 X2 ■ (C2) 

The shock front travel speed derived from the upstream and 
downstream sides should be the same. Thus, from expres- 
sions <|Cljl and l|C20 , we derive the discontinuity in constant 
k, namely 

(fe 2 Ai) 1/2 = xi/xa . (C3) 

This discontinuity in k originates physically from the dis- 
continuity in temperature or entropy across the shock front, 
as discussed in Section 2.3 of the main text. We also list flow 
velocities in the shock reference framework here 

iti — ui — Us — k^^t"^ 1 (vi — nxi) (C4) 

and 

U2 = U2 — Its = k)/ 2 t n ~^(v2 ~ UX2) ■ (C5) 

Referring to the sound speed definition J8J, the pre- 
shock (upstream) sound speed is c\ — r y 1 ^ 2 kl^ 2 t n ~ 1 a\ y ~ 1 ^ 2 . 
Then from equations 1171 and 1C4L the pre-shock (up- 
stream) Mach number is 

r. _ Ui _ kl^t"' 1 ^! - nxj) v x - nx\ 

" Cl " fc^V-V/2 a C7-l)/ 2 - 7 l/2 a (7-l)/2 ' W 

where M is the Mach number in the reference framework of 
the shock front. With this derived upstream Mach number, 



